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ABSTRACT 

In direct sequence-code division multiple accesses systems (DS-CDMA), the signals are transmitted over 
multipath channels that introduce fading. Multipath fading along with multiple access interference and inter - 
symbol interference degrades the system performance. Further, simultaneous presence of multipath fading and 
shadowing leads to worsening of wireless channels. Moreover, experimental results have confirmed the 
presence of impulsive noise in wireless mobile communication channels. Hence, this paper presents a 
technique for multiuser detection in DS-CDMA systems over shadowed fading channels in presence of 
impulsive noise. Approximate expression for average probability of error of an M-decorrelator is derived, for 
the demodulation of binary phase shift keying (BPSK) signals over shadowed fading channels, by modeling 
the channel with generalized K (GK) distribution. A new M-estimator is proposed for robustifying the detector 
and its performance is also studied and analyzed by evaluating probability of error with the derived 
expression. Simulation results reveal that the proposed multiuser detector performs better in fading and 
shadowing with heavy-tailed impulsive noise. 
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I. INTRODUCTION 

Recent research has explored the potential benefits of 
evolutionary optimization algorithms and their 
application to multiuser detection (MUD) for direct- 
sequence code division multiple access (DS-CDMA) 
systems [l]-[4]. An adaptive robust MUD technique 
for CDMA by implementing Ruber's M-estimator 
using genetic algorithm (GA) is presented in [3] and 
shown that it is robust against heavy-tailed 
impulsive noise. Recently, particle swarm 
optimization (PSO) algorithm has been applied for 
the MUD [4] to detect received data bit by optimizing 
an objective function. 

The generalized-K (GK) fading channel has been 
received considerable attention as it can provide a 
good fit to different fading environments such as 
Nakagami-m and Rayleigh-Lognormal [5]. 
Experimental results have confirmed the presence of 



heavy-tailed impulsive noise in outdoor mobile 
communication channels, in radar and sonar systems 
and in indoor wireless communication channels [6] . 
Hence, this paper presents the implementation and 
performance analysis of proposed based M-estimator 
[7], [10], which performs well in the heavy-tailed 
impulsive noise, using PSO algorithm. 

II. SYSTEM MODEL 

Consider an L-user synchronous CDMA system, 
where each user transmits information by 
modulating a PN sequence over a single-path GK 
fading channel. The received signal over one symbol 
duration can be modeled as [8] 

MO ^ 



L M-1 

rit) = ^\T. S bi{i)aiit)Q^'^''^Siit-iT^-Ti)\ + n{t) 
1=1 i^O 



(1) 



where 9^ {•} denotes the real part, M is the number of 

data symbols per user in the data frame of interest, Ts 
is the symbol interval, a/(0 is the time-varying 
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fading gain of the 1*^ user's channel, ^(f) is the time- 
varying phase of the 1*^ user's channel, Z?/ (/) is the i*^ 

bit of the 1th user, ^^(Ois the normalized signaling 
waveform of the 1*^ user and n(t) is assumed as a 
zero-mean complex two-term non-Gaussian noise 
[9]. 

For synchronous case (i.e., = r2 = ... = = 0), 

assuming that the fading process for each user varies 
at a slower rate that the magnitude and phase can 
taken to be constant over the duration of a bit, the 
received signal can be expressed in matrix notation 
as [9] 

r(i) = AOii) + w{i) (2) 
where r(/) □ [r^ (/),... ,r^(/)] , 

T 

□ [ivi (/),.. and 

^(OD ^[^i(0^i(0,...,^L(0gL(0f . 

Here, >v^(/) is a sequence of independent and 

identically distributed (i.i.d.) complex random 
variables whose in-phase and quadrature 
components are independent non-Gaussian random 

variables, gi (i) is the 1*^ channel fading coefficient 
and AD [^^,^2,...,^^] witha^ □ [a|,a2,...,a|y^]. 

It is assumed that the signal of each user arrives at the 
receiver through an independent, single-path fading 
channel. For the shadowed fading channels, a^ii) are 

i.i.d. random variables with GK distribution given by 
[12] 
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where, m is the Nakagami fading parameter that 
determines the severity of the fading, //represents 

the shadowing levels, Q.^ is the average SNR in a 

shadowed fading channel, K^{-) is the modified 

Bessel function and r(-) is the Gamma function [9] 



In M-estimates, unknown parameters , ^2 ' •••^L 
solved by minimizing a sum of function, p{-) of the 
residuals [9] 
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where p{') represents a specific penalty function that 

is symmetric positive-definite with a unique 
minimum at zero, and is chosen to be less increasing 

than square, r^{i) and (/) are the n^^ ^^vd 
elements of the vectors r(/) and 9{i) respectively, 
[4]^/ is the nl^h element of the matrix A , and 

denotes imaginary part. This paper considers 

the modified-Hampel based M-estimator to 
implement an M-decorrelating detector with penalty 
function [7] 
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where a, b and d are constants that depend on the 
robustness of the estimator [7]. 

III. PSO BASED M-DECORRELATING DETECTOR 

PSO is a swarm intelligence method for global 
optimization modeled after the social behavior of 
bird flocking and fish schooling [4]. In PSO 
algorithm the solution search is conducted using a 
population of individual particles, where each 
particle represents a candidate solution to the 
optimization problem (3). Each particle keeps track 
of the position of its individual best solution (called 

1 , s best r best best n j.i 11111 

as pbest), =VP^^ ,...,p^^ J and the overall global 

1 , 1,' /11J 1 i\ best r best best 

best solution (called as gbest), g =[g ,...,g J 

among p bests of all the particles in the population 
represented as = , p^^^ ] , where p ^ is the d 



th 
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particle in the it^^ iteration, d = 1, 2, . . ., Np, it = 1, 2,. . 
Nit, Np is the number of particles, and Nit is the 
maximum number of iterations. Corresponding to 
each position, the particle velocity is 

v^^ = [v^ v^^ ] [2], [4]. The steps involved in the 

PSO based M-decorrelating detector's 
implementation are [2], [3], [4]: 

Step 1: Compute the decorrelating detector output, 

^ =R-^£^r. Here, R (□ AJn) is the 

normalized cross-correlation matrix of signature 
waveforms of all users. 
Step 2: Initialization: The output of decorrelating 

detector is taken as input first particle = 9j . 

Step 3: Fitness evaluation: The objective function (3) 
is used to find the fitness vector by substituting 

residuals. Local best position p^^^^is recorded 

by looking at the history of each particle and the 

particle with lowest fitness is taken as g of 

the population. 
Step 4: Update the inertia weight by using the 

decrement function w^^ - Pw^ ^ , where y^<lis 

the decrement constant. 
Step 5: Update the particle velocity by using the 
relations 

privatethinking of the particle 

it+l it 

(6) 



it+l ... 



( best it \ 

(best it \ 



■V- 



social thinking of the particle 



it + l it it + l /ry\ 



where c^ and c^ are the acceleration constants 

representing the weighting of the stochastic 
acceleration terms to pull the particle to pbest and 

gbest. ju^ and lu^ are random numbers that are 

uniformly distributed between 0 and 1. Particle 
position is updated according to (6). Particle velocity 
is limited by the maximum 

1 . , max r max max 

velocity V = [v^ .--.v^ J. 



Step 6: The individual best particle position is 



updated by following rule: 



, -I best it 

then = p^ 



best 



Step 7: g is the global best particle position among 
all the individual best particle positions p^^ at 

the it*h iteration such that F {^g^^^ ) - ^ (p^J ) ' 

Step 8: The above steps are repeated until the 
maximum number of iterations has been 
reached. 

The computed g^^^^ value is used to evaluate average 

probability of error of BPSK demodulator, over GK 
fading channel, by an approximate expression given 
by 
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where m is the Nakagami fading parameter that 
determines the severity of the fading, jU represents 

the shadowing levels, is the average SNR in a 



shadowed fading channel F = 



\mju 
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and 



mju 



and ^ (•) is the Whittaker function [13]. 



o 



TABLE I. PSO PARAMETERS USED FOR SIMULATION 



Parameter 


^alue 


Number of particles, Np 


20 


Maximum number of iterations. Nit 


100 


Acceleration constants and 


2 


Maximum velocity of the particles 
V^"^^ for all users 


2 


Initial inertial weight, w 


1 


Decrement constant, J3 


0.99 
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IV. SIMULATION RESULTS 

In this section, the performance of M-decorrelating 
detector is presented by computing (7) for for 
different values of fading parameters and shadowing 
levels with least-squares, Huber, Hampel and 
modified-Hampel penalty functions. The PSO 
parameters used for simulation are presented in 

Table L The value of g^^^^ is computed for specified 

penalty functions and is used to compute the average 
probability of error (7) of BPSK signals. 

In Fig. 1 and Fig. 2, the average probability of error 
versus the signal-to-noise ratio (SNR) corresponding 
to the user 1 under perfect power control of a 
synchronous DS-CDMA system with six users (L = 6) 
and processing gain, N = 31 is plotted for moderate 
impulsiveness {s= 0.01) of noise, m = 1, |a = 1. 
Similarly, the average probability of error is plotted. 

In Fig. 3 and Fig. 4 with highly impulsive noise ( e = 
0.1) and m = 1, \x = l. These simulation results show 
that the decorrelating detector with proposed 
estimator based detector performs well even in 
highly impulsive noise when compared to least- 
squares, Huber and Hampel estimator based 
detectors. 

V. CONCLUDING REMARKS 

Multiuser detection technique for DS-CDMA systems 
over GK fading channels using PSO was presented in 
this paper. An objective function, which is a sum of 
less rapidly increasing function of residuals, was 
used to obtain optimum estimates. An M- 
decorrelator is implemented with different influence 
functions. Simulation results shows that the 
proposed M-decorrelator performs better when 
compared to least-squares, Huber and Hampel 
estimator based detectors. 
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Fig.l Average probability of error versus SNR for 
user 1 for linear multiuser detector, minimax detector 
with Huber, Hampel and proposed M-estimator in 
synchronous CDMA channel with impulse noise, N = 
31, £ = 0.01, m = 1, |a = 1. 



jj epsilon-O.Ol, kappa-100. L-6 
10 £ 1 \ 1 1 1 1 1 r 




SNRuulD 

Fig.2 Average probability of error versus SNR for 
user 1 for linear multiuser detector, minimax detector 
with Huber, Hampel and proposed M-estimator in 
synchronous CDMA channel with impulse noise, N = 
31, £ = 0.01, m = 2, \i = 2. 
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Fig.3 Average probability of error versus SNR for 
user 1 for linear multiuser detector, minimax detector 
with Huber, Hampel and proposed M-estimator in 
synchronous CDMA channel with impulse noise, N = 
31, 6 = 0.1, m = 1, |a = 1. 
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Fig.4 Average probability of error versus SNR for 
user 1 for linear multiuser detector, minimax detector 
with Huber, Hampel and proposed M-estimator in 
synchronous CDMA channel with impulse noise, N = 
31, £ = 0.1, m = 2, \i = 2. 
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